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ABSTRACT. 

We present a self-consistent nonparametric model of the cosmic velocity field 
based on the spatial distribution of IRAS galaxies in the recently completed 
all— sky PSCz redshift survey. The dense sampling of PSCz galaxies allows us 
to infer peculiar velocities field up to large distances with unprecedented high 
resolution. 

The most streaking feature of the PSCz model velocity field is a coherent 
large— scale streaming motion along the Perseus Pisces, Local Supercluster, 
Great Attractor and Shapley Concentration baseline, with no evidence for 
a backinfall into the Great Attractor region. Instead, material behind and 
around the Great Attractor in inferred to be streaming towards the Shapley 
Concentration. 

A likelihood analysis that uses the information available on bulk velocities, 
cosmological dipoles and local shear has been performed to measure (3. We 
have obtained (3 = 0.6^q' 22 (1-cr), in agreement with other recent determina- 
tions. 



1 Introduction 

The Gravitational Instability pic- 
ture and the Linear Biasing hypoth- 
esis provide us with a theoretical 
framework in which galaxy peculiar 
velocities are related to the inhomo- 
geneities of the underlying mass den- 
sity field, 8 m . On large scales, where 
Linear Theory applies, the relation 
is remarkably simple and one can 
compute peculiar velocities directly 
from the galaxy distribution. The re- 
sulting velocity field is fully speci- 
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fied by the (3 = " parameter, 
where b is the bias parameter that 
linearly relates the inhomogeneities 
in the galaxy distribution to the un- 
derlying mass density field. Compar- 
ing the model velocities with the ob- 
served ones allows one to check the 
plausibility of the adopted theoreti- 
cal scenario and to constrain the /3 
parameter. Here we present a new 
model for the velocity field predicted 
from the PSCz survey that we com- 
pare with existing peculiar velocity 
catalogues. 



2 The PSCz Galaxy Survey 

Large sky coverage is the basic 
requirement to model peculiar ve- 
locities from a redshift survey. The 
dataset used in this work is the re- 
cently completed PSCz redshift sur- 
vey described in detail by Saunders 
et al. (1998) and briefly summarized 
here. The main catalogue contains 
some 15,500 IRAS PSC galaxies with 
60 /im flux, /60, or greater. To avoid 
cirrus contamination only PSC ob- 
jects with /xoo < 4/eo were se- 
lected. Stars were excluded by re- 
quiring that feo > 0.5/25- For our 
purposes, one of the most impor- 
tant properties of the PSCz cata- 
logue is its large sky coverage. The 
only excluded regions are two thin 
strips in ecliptic longitude that were 
not observed by the IRAS satellite, 
the Magellanic clouds, and the area 
in the galactic plane where the B— 
band extinction exceeds 2 magni- 
tudes. Overall, the PSCz catalogue 
covers ~ 84% of the sky. It is worth 
stressing that the survey is deeper 
and the galaxy sampling is denser 
that in any previous all-sky cata- 
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Figure 1. Distribution of IRAS PSC galaxies vs. redshift. The upper and 
the lower, shaded histograms represent the distribution of PSCz and 1.2Jy. 
galaxies, respectively. The curves show the expected counts as a function of 
distance estimated from the selection functions. The heavy line at the bottom 
shows the predicted distance distribution of Abell/ACO clusters. The labels 
give the total number of objects in each sample. 



logue. This point is clearly shown 
in Figure 1 in which the histogram 
PSCz galaxy counts as a function of 
redshift is compared to the 1.2 Jy 
one (Fisher et al. 1995) and to the 
expected distribution inferred from 
the luminosity function (continuous 
lines) . 

3 Building Velocity Models 

Systematic errors are a major con- 
cern when modeling, measuring and 
comparing peculiar velocities. In the 
modeling phase we keep systematics 
under control by using three inde- 
pendent techniques that predict ve- 
locities from the galaxy distribution 
in redshift space: 

1) A particle-based iterative 
method similar to that introduced 
by Yahil et al (1991), which pre- 
dict the peculiar velocity of the 
PSCz galaxies at their real space 



positions. 

2) A similar iterative algorithm in 
which, however, densities and ve- 
locities are computed onto a reg- 
ular grid. 

3) The Spherical Harmonics ex- 
pansion technique introduced by 
Nusscr and Davis (1994) in which 
the velocity field is predicted at 
any points in the redshift space. 

All these methods assume linear or 
quasi-linear theory and therefore re- 
quire some degree of smoothing. 

Random and systematic errors for 
the three model velocity fields have 
been evaluated using a suite of mock 
PSCz catalogues extracted from the 
N-body simulations by Cole et al. 
(1997). 

4 A Cosmographic Tour 

The depth and sampling frequency 
of the PSCz datasct allow a reliable 
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Figure 2. PSCz density and velocity fields within 120 h —1 Mpc on a slice 
through the Supergalactic plane. Continuous lines represent positive isoden- 
sity contours drawn with a spacing of A<5 m = 0.5. The thick line sets the 
overdensity level. Negative contours are drawn with a dashed line. The am- 
plitude of the velocity vectors, obtained for (3 = 1, is on an arbitrary scale. 



map of the density and velocity fields 
to be constructed within a sphere 
around us of 120 /i -1 Mpc. Figure 2 
shows the PSCz model density and 
velocity fields smoothed with a 6 ft -1 
Mpc Gaussian filter in a slice through 
the Supergalactic plane. 

The Local Supercluster, cluster at 
(SGX,SGY) = (-2.5,11.5), is con- 
nected to the 
prominent Hydra-Ccntaurus super- 
cluster at (SGX,SGY) = (-35,20). 
Together with 
the Pavo— Indus-Telescopium super- 
cluster [(SGX, SGY) = (-40,-15)], 
the latter makes up the well known 
Great Attractor. The Coma cluster 
appears as a peak at (SGX, SGY) = 
(0,75). The Perseus-Pisces super- 
cluster is at (SGX, SGY) = (45, -20) 



and its northern extension is visible 
at [(SGX, SGY) = (45,20)]. The Ce- 
tus Wall may be seen as an elongated 
structure around (SGX, SGY) = 
(15, -50). Finally, the Shaplcy Con- 
centration starts appearing towards 
the edges at [(SGX, SGY) = 
(-90,60)]. The Sculptor void is 
the largest underdense region at 
[(SGX, SGY) = (-20,-45)], almost 
matched in size by the void centered 
at [(SGX, SGY) = (70, 50)]. 

The local velocity field implied by 
the density field is dominated by 
the large infall patterns towards the 
Great Attractor, Perseus-Pisces and 
Coma. A striking feature is the large- 
scale coherence of the velocity field, 
apparent as a long ridge between 
Cctus and Perseus-Pisces and as a 
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Figure 3. Model vs. observed bulk velocity. Amplitude of the bulk flow vectors 
(lower left-hand) and of its SGY -component (upper left-hand) are shown. 
The various symbols represent different experimental determinations. The 
dashed lines indicate the 1-a error strip from PSCz prediction. The vectors' 
directions, measured at 50 h~ x Mpc, are displayed in the right-hand plot. The 
filled circle represents the prediction from PSCz. 



large-scale flow along the Perseus- 
Pisces (North), Virgo, Great Attrac- 
tor, Shapley Concentration baseline. 
Note that a firm prediction of the 
PSCz data is the lack of prominent 
back-infall onto the Great Attractor 
region. 

5 The Bulk Flow 

The bulk velocity represents one 
of the simplest low-order statistics 
that, in principle, can be estimated 
observationally and which has a the- 
oretical counterpart. Measuring and 
modeling the bulk flow within a 
given volume, however, is prone to 
random and systematic uncertain- 
ties. Both errors have been estimated 
using the mock PSCz catalogues. In 
Figure 3 we compare our theoret- 
ical predictions to various observa- 
tional measurements. On the left- 
hand side we show the amplitude of 



the bulk velocity (lower plot) and 
of its SGY component (upper plot). 
The strip delimitated by the dashed 
lines represents the model predic- 
tion and its 1-a uncertainties. Filled 
triangles are taken from Dekel et 
al. (1998) and represent the bulk 
flow from the Mark III catalogue 
(Willick et al. 1997), measured using 
a POTENT-smoothing technique. A 
similar method has been applied 
by Eldar et al. (1998) on the SFI 
catalogue (Giovanclli et al. 1997). 
Their result is displayed with filled 
square. The open triangle shows the 
bulk flow obtained by considering the 
sample of SNc la (Riess et al. 1995). 
The open square displays the result 
by Lauer and Postman (1994). The 
right-hand side of the figure shows 
the direction in the sky of the vari- 
ous bulk flows. The filled circle rep- 
resents the PSCz prediction. The di- 
rection of the CMB dipolc is plotted 
for reference and is represented by an 
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asterisk in the centre of the figure. 

The bulk flow predicted using 
PSCz agrees well, both in amplitude 
and directions, with all the experi- 
mental determinations but, perhaps 
not surprisingly, disagrees with the 
result obtained by Lauer and Post- 
man. Comparing predicted and mea- 
sured bulk flow allows one to mea- 
sure (3. As indicated in the figure, we 
find P ~ 0.75 with, however, a large 
error of ~ 50 % 



6 Estimating (3 

Predicted peculiar velocities can 
be compared to the observed ones to 
estimate j3. This exercise, however, 
is potentially prone to systematic bi- 
ases that need to be corrected for. We 
have performed a likelihood analysis, 
similar to that introduced by Strauss 
et al. (1992), in which such a com- 
parison can be performed that auto- 
matically accounts for random and 
systematic errors in a statistical fash- 
ion. The observational quantities we 
consider are: the velocity of the Lo- 
cal Group measured from the CMB 
dipole (vlq = 625 ±25 km s _1 from 
the 4 years COBE data of Linewaver 
et al. (1996)), the bulk flow measured 
from the Mark III catalogue by Dekel 
et al. (1998) and the evidence of a 
small "shear" (|v s | < 200 km s _1 ) 
indicating that the local flow is re- 
markably cold (see van de Weygaert 
this volume). The theoretical quanti- 
ties, predicted from the distribution 
of PSCz galaxies, are: the accelera- 
tion of at the position of the Local 
Group and the cumulative bulk flow 
in spheres of increasing radius. Un- 
der the hypothesis that all the quan- 
tities mentioned above are Gaussian 
random fields and assuming a prior 
model for the underlying cosmology 
it is possible to obtain an analytic 
expression for their joint probability 
distribution. Maximizing the distri- 
bution with respect to the free cos- 
mological parameters allows one to 
estimate (3. We have obtained that 
/3 = 0.6±°' (1-cr). Along with the 
constraints that this analysis sets on 
the other parameters, we conclude 
that a universe with high f2 m is less 
favourite by observations. 

Other techniques of constraining 
(3 involve point-by-point velocity 
comparisons. Because of its dense 



sampling and low shot noise errors 
the PSCz velocities can be reliably 
modeled out to large distances and 
then compared with objects capa- 
ble of tracing the large scale dynam- 
ics, such as clusters of galaxies. A 
velocity— velocity comparison of this 
type has been recently performed us- 
ing the PSCz model velocity field 
and the SMAC catalogue of clusters' 
velocities and is described by Hudson 
et al. in this volume. 
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